Spin-glass state and long-range magnetic order in Pb(Fei/2Nbi/2)03 
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We have investigated the magnetic ground-state of the multiferroic relaxor ferroelectric 
' Pb(Fei/2Nbi/2)03 with /iSR spectroscopy and neutron scattering. We find that a transition to 

. a partially disordered phase occurs below T — 20 K that coexists with long-range antiferromagnetic 

' ordering. The disordered phase resembles a spin-glass. No clustering of magnetic ions could be 

I evidenced by ^SR so that the coexistence appears homogeneous in the sample. 
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PACS numbers: 75.50.Lk Spin glasses and other random magnets; 77.80.-e Ferroelectrieity and antiferroelec- 
tricity; 61.05.F Neutron diffraction and scattering; 76.75.+i Muon spin rotation and relaxation 



I. INTRODUCTION 



c/2 I Since the discovery of a large magneto-electric effect in RMn203EI and RMn205H (R=Rare Earth), there is a revival 
of interest in the study of multiferroic materials where the coupling between the electric polarization and the magnetic 
order is strong. Although coexistence of ferroelectrieity and magnetic long-range order is, not very common, quite a 
number of new multifcrroics materials have been found recently, like Ni3V208,El MnW04,a or LiCuV04j3 to cite only 
a few. Multifcrroics that have coupled ferroelectric and magnetic order parameters are also promising materials for 
applications, as it is possible to control the ferroelectric polarization by a magnetic field and vice-versa.El Therefore, 
there is an extensive search for materials that show pronounced magneto-electric effect at ambient conditions. B-site 
I , disordered complex perovskites PbBfpB"_^03 might be attractive candidates as they can adopt various ions which 
''O ' strongly influences the temperature of magnetic and ferroelectric ordering. 
^ I Lead iron niobate, Pb(Fei/2Nbi/2)03 (PFN ), is one of the well known complex perovskites that exhibits multiferroic. 
^ ■ properties Ji PFN is rhombohcdral at room temperature and cubic above the ferroelectric Curie temperature of 383 K.eI 
I i] Pb(Fei/2Nbi/2)03 undergoes a diffuse phase transition without pronounced frequency dispcrsion.EI In the cubic phase, 
the symmetry of the chemical cell is Pm3m. There is a cubic-tetragonal phase transition at T = 383 K followed by 
a first-order transition at 354 K. The syruHietry of the chemical structure of the low-tcmpcraturc phase of PFN was 
^ , the subject of. some debate. Lampis et alS3 reported that the structure is monoclinic with space group Cm, whereas 
Ivanov et al.ll using powder neutron diffraction data found that the chemical structure of Pb(Fei/2Nbi/2)03 has i?3c 
QQ , symmetry at T = 10 K. 

^-H ■ The magnetic properties of PFN are also notxpmpletely understood. PFN undergoes a paramagnetic to antiferro- 
' I magnetic phase transition at 140-160 K.tj The magnetic structure was determined by single crystal neutron 

. diffraction and corre^onds to a simple G-type antiferromagnetic arrangement of the magnetic moments at all tem- 



perature below Tat.Ej Magneto-electric effect in PFN at low temperatures has been demonstrated long ag oy and, 
. recently, iLwas shown that the DC magnetic susceptibility has an anomaly at the paraelectric-ferroelectric phase 
' transition.li3 The dielectric permeability also shows an anomaly at the Neel temperatureE^I, which suggests that a 
. ^ , biquadratic term in the Landau free energy of the form P^M^ is responsible for the coupling between the ferroelectric 
' polarization P and the (staggered) magnetization Af . Additional magnetk-|anomalies were reported in the lilprature 
JH , for Pb(Fei/2Nbi/2)03 below the magnetic phase transition. Both NMR liZI and susceptibility measurementaia have 
shown that an additional magnetic phase transition occurs in PFN around Tg ~20 K that was interpreted as the 
formation of a spin-glass phase by Kumar ct al.tj 

Here we report on muon-spin rotation and neutron scattering results of the magnetic properties of PFN below 
the Neel temperature. We show that the magnetic ground-state of Pb(Fei/2Nbi/2)03 is a spin-glass-likc state that 
coexists with long-range antiferromagnetic order below Tg c:^ 20 K. 



II. EXPERIMENTAL 



The Pb(Fei/2Nbi/2)03 single crystals used in the muon spin rotation and neutron scattering experiments were grown 
by spontaneous crystallization from the melt following the routine described in Ref. |l^. The muon spin relaxation 
measurements were performed using the GPS instrument at the Paul Scherrer Institut (Villigen, Switzerland) on a 
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Pb(Fei/2Nbi/2)03 single crystal between 160 K and 5 K. The data were recorded using the zero-field method that 
allows to determine both the static and dynamics in disordered spins sjistems. The neutron scattering experiments 
were performed with the cold neutron three-axis spectrometer TASPE3 at SINQ operated in diffraction mode at 
kf =1.97 A~^. For this experiment a single crystal of Pb(Fei/2Nbi/2)03 was oriented with the [1,1,0] and [0,0,1] 
crystallographic directions in the scattering plane. In that scattering geometry, magnetic Bragg peaks with indexes 
{h/2, h/2, 1/2) in the rhombohcdral setting can be accessed. The use of a three-axis spectrometer was justified by the 
need to improve the ratio intensity/background in order to search for possible short-range order at low temperatures. 



III. RESULTS 



A. ^SR 

In a /iSR experiment, the time evolution of the muon spin polarization is monitored by recording the asymmetric 
emission of positrons produced by the weak decay of the muon (muon lifetime ~ 2.2 fis). The time histogram of the 
collected positrons is given by 

N{t) = No exp{-t/T) [1 + AG^ (t)] +B (1) 

where A is the initial muon asymmetry parameter, A'o is a normalization constant, and i? is a time-independent 
background. The function Gz{t) reflects the normalized muon-spin auto-correlation function 

(5(^:5(0)) 

52(0) 

where S is the spin of the muon. Hence AGz{t)., often called the jjSK signal, reflects the time evolution of the muon 
polarization. 

Figure |l] shows the ^SR signal in Pb(Fei/2Nbi/2)03 upon passing the Neel temperature. In the paramagnetic phase 
the signal is consistent with a weak depolarization of the muon signal solely due to nuclear dipole moments. Below the 
Neel temperature, the muon signal changes rapidly and the initial asymmetry A is strongly reduced. This indicates 
that there is a large distribution of internal fields seen by the which is a consequence of the random distribution 
of the Fe and Nb ions in the chemical lattice. Therefore, the /iSR signal is depolarised and the residual asymmetry 
is close to 1/3 of the value observed in the paramagnetic state. Hence, in the long-range antiferromagnetic phase, 
the muon signal reflects only muons having their initial polarization along the direction of the internal field at the 
stopping site. 

Below Tjv, the time evolution of the remaining /iSR signal is best fitted assuming the function 

G,(t) = iexp[-(Ai)^]. (3) 

Close to , the power exponent is /3 ~ 1 , reflecting the presence of fluctuating internal magnetic fields with a unique 
correlation time Tc- In this range of temperatures, the depolarization rate A is proportional to the second moment of 
the magnetic field distribution {B"^) and to the correlation time of the internal field and is given by 

A = 7^(SVc, (4) 

where 7^ is the muon gyromagnetic ratio. The temperature dependence of the depolarization rate A is shown in 
Fig. H It can be observed that there is a critical- like divergence close to Tg = 20 K. This is is a clear evidence that 
Pb(Fei/2Nbi/2)03 undergoes a second magnetic phase transition at the temperature Tg. 

Figure ^ depicts the temperature evolution of the /3 exponent. As already mentioned, at temperatures just below 
T/v, the muon depolarization is described by an exponential function instead of a stretched exponential, i.e., with 
the parameter (3 ^1. Upon further cooling below T^r, (3 continuously decreases and approaches the value of 1/3 at 
Tg ~ 20 K. A decrease of the exponent (3 toward the value, of 1/3 has been predicted theoretically and later observed 
for the dense spin-glass systems AgMn and AuFe alloys.Eil It is tempting to attribute this change of (3 to the change 
in the distribution of the correlation times Tc of the internal magnetic fields at the muon stopping sites as a spin-glass 
transition is approached from above. However, for an usual paramagnetic to spin-glass transition, one observes that 
the residual fluctuations lead to a constant value of (3 below Tg while A decreases, reflecting the decrease of the 
magnitude of the fluctuating field, as the temperature is lowered deeper into the spin-glass phase. However, Fig. ^ 
also shows that (3 grows as the temperature is lowered below Tg to recover the value ~ 1 at T = 2 K. The solution 
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of this apparent paradox is that the slow spin fluctuations, arising from the spin-glass state and which persist in a 
relatively broad temperature range below T^, occur on top of antiferromagnetic fluctuations due to the long-range 
order. It is then natural to expect that the strong reduction of the critical spin-glass fluctuations in time when the 
temperature is decreased below Tg leads to the recovery of the weak non-critical antiferromagnetic fluctuations probed 
by the /xSR signal above Tg. We note that as the static component of the /iSR signal is already lost above Tg, no 
indication ota change (i.e., increase) of the damping of this component (as evidenced for example in other reentrant 
spin-glasseg23) could be determined at the spin-glass transition by /xSR. 



B. Neutron scattering 

Figure ^ shows a typical neutron elastic scan along the (q, q, —2q) direction at T = 5.2 K in Pb(Fei/2Nbi/2)03 . The 
spectrum consists of a sharp Bragg peak centered around Qjv =r(f /2, 1/2, 1/2) that corresponds to the antiferromag- 
netic long-range order detected by neutron diffraction previouslytj. A diffuse component is also present in the neutron 
spectrum that is weaker in intensity and has a width in q that is much broader than the experimental resolutiop. The 
magnetic origin of the diffuse component was confirmed using the neutron polarization analysis device MuPADlj that 
allows a complete separation of magnetic and nuclear scattering. Namely, using a neutron polarized beam for which 
the polarization is aligned along the scattering vector, magnetic scattering reverses the neutron polarization vector, 
whereas scattering of nuclear origin leaves the polarization of the neutron beam unchanged. Magnetic scattering 
therefore appears in the spin-flip channel, as shown to be the case for the diffuse component in the inset of Fig. ^. 

We investigated the distribution of the diffuse scattering around the Qjv antiferromagnetic Bragg peak and found 
that it is essentially isotropically distributed in the (1,1,0)/(0,0,1) scattering plane accessible in the present measure- 
ments (see Fig. |^) . We performed an energy scan through the diffuse component and found that it is resolution-limited 
within the energy window {AE=0.3 meV) of the spectrometer. Therefore, the spin fluctuations associated with the 
diffuse scattering are static in the the time-window of our experiment. 

Within the quasi-statiCpapproximation, the energy-integrated neutron cross section is proportional to the generalized 
static susceptibility x(Q)Ej 

/(Q, T) ^Tj2iSc.p - E FdiQ)Fd> {Q)xii(.Q) 

Q/3 dd' ^r^-j 

X exp (iQ{rd - r'^)) . 



The term (^a/3 — QqQ/s) accounts for the fact that only fluctuations perpendicular to the scattering vector Q contribute 
to the diffuse intensity. The index d labels the spins at position and Fd{Q) is the magnetic form factor. We modeled 
the diffuse scattering measured along {q,q, —2q) around the Qjv magnetic Bragg peak by a Lorentzian function 

^(Q:^)- 2^io^^n 12 (6) 

where A{T) is the scattered intensity and k = l/'C is the inverse of the correlation length. The temperature dependence 
of the diffuse scattering intensity is shown in Fig. Above T ~120 K the diffuse scattering is too weak to be separated 
from the background scattering. Below that temperature A{T) increases monotonically with cooling down to T ~2 K. 
A flt to the neutron data with Eq. (^ yielded at all temperatures, for which the diffuse scattering could be observed, 
K 0.015 (rlu) that corresponds to ^ ^ 17 A. 



IV. CONCLUSION AND DISCUSSION 



Both /iSR spectroscopy and neutron diffraction have shown that a transition to a spin-glass state that coexists with 
long-range antiferromagnetic order occurs in the rclaxor pcrovskitc Pb(Fei/2Nbi/2)03 a-t Tg ~ 20 K. Figure]^ shows 
that the asymmetry of the muon signal is rapidly lost below T^ and does not recover when passing the spin-glass 
transition. This indicates that the magnetic-order is homogeneous and that it persists upon cooling the sample below 
Tg. This is confirmed by the presence of the resolution-limited magnetic Bragg peak at the lowest temperature showing 
that long-range antiferromagnetic order of the Fe'^^ spins is still present below Tg. The long-range antiferromagnetic 
order coexists with the disordered spin-glass-like phase. The signature of the spin-glass order appears in the neutron 
spectrum in the form of diffuse scattering. 
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There arc numerous examples of coexistence between the ferromagnetic and spin-glass states of matteii^^. For exam- 
ple, ferromagnetic and spin-glass orders are known to coexist in AuFe alloys at sufficiently low temperaturesc^l. Antifet.- 
romagnetic and spin-glass orders are also known to coexist at sufffcicntly low temperatures, e.g., in Feo.55Mgo.45Cl2E3. 
The coexistence of collinear magnetic long-range order, e.g., ferromagnetic or antiferromagnetic, and spin-glass order 
is well understood theoretically at the level of the Sherrington-Kirkpatrick (SK)ca infinite-range model for classical 
(Ising- or Heisenberg) spins. In the SK model, this coexistence is the rule below some critical temperature and down 
to the lowest temperatures as soon as the magnitude | J| of the mean value of the random exchange counIing^|L,hiitwcen 
any two spins is of the order of or larger than the standard deviation A J of its Gaussian distributionEj'EHlLilH 

Classical spitunodels with short-range but random two-body spin interactions (the Edwards- Anderson (EA) model 
of spin glassea^j) are not as well understood as the SK model. Monte Carlo simulations suggest that the lower 
critical dimension for the cpset of spin-glass order when the disorder is dominant, | J|/AJ -C 1, is strictly below three 
dimensions for Ising spins,Eil a prerequisite for the coexistence of collinear and spin-glass order when | J|/AJ > 1x3 
Fluctuation effects are stronger for Heisenberg spins than they are for Ising spins, as fluctuations in the directions 
transverse to the direction taken by the collinear order parameter in spin space become available in the former case. 
Correspondingly, state of the art Monte Carlo simulations give "only" the uppet bound of three dimensions for the 
spin-glass lower critical dimension when the disorder is dominant, \ J\/AJ ^ lE3 In other words, it cannot be ruled 
out, within the accuracy of available Monte Carlo simulations, that the spin-glass phase is quasi-long-range ordered 
instead of long-range ordered when the classical spins are Heisenberg-like in the three-dimensional EA model. 

Given the uncertainty of the spin-glass lower-critical dimension for Heisenberg EA spin-glass models, the mere 
observation that antiferromagnetic and spin-glass orders coexists in Pb(Fei/2Nbi/2)03 is not sufficient to establish 
whether the relevant magnetic degrees of freedom are Ising- or Heisenberg-like. Although exchange interactions 
remain to be determined in Pb(Fei/2Nbi/2)03 j the relatively high Neel temperature of Pb(Fei/2Nbi/2)03 suggests 
that exchange interactions dominate dipolar and single-ion anisotropics in this compound. If this is the case, then 
Pb(Fei/2Nbi/2)03 is a good candidate to study the coexistence of long-range antiferromagnetic and spin-glass orders 
for the Heisenberg spins of a three-dimensional antiferromagnetic EA model. 
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Figure 1: (Color online:) (a) Time evolution of the muon-spin polarization reported for characteristic temperatures. 

Note the weak depolarization in the paramagnetic phase (160 K) and the loss of the signal amplitude below T/v, 
which retains the value 1/3 {t = 0) at low temperatures, (b) Examples of polarization signals recorded close to (20 
K) and well below (5 K) Tg. The lines represent fits performed with Eq. (|[) (data reported on (b) are strongly 

binned for clarity). 
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Figure 2: Temperature dependence of the ^SR depolarization rate A below 80 K. 
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Figure 3: Temperature dependence of the parameter /3, showing the occurrence of a distribution of relaxation times 

in PFN (see text). 
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Figure 4: (Color online:) Transverse scans through the Qn = (0.5, 0.5, 0.5) magnetic Bragg peak. The broad 
component in the spectrum was fitted by a Lorentzian and is emphasized by a bold line. The Bragg peak is a fit to a 
Gaussian. The inset shows the intensity recorded in the neutron spin-flip channel for incident polarisation along the 

scattering vector. 
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Figure 5: (Color online:) Distribution of the elastic diffuse scattering measured in Pb(Fei/2Nbi/2)03 at T = 2 K. 
For clarity the intense Bragg peak at Qtv = (0.5, 0.5, 0.5) was removed and the background subtracted from the data. 
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Figure 6: Temperature dependence of the integrated intensity of the diffuse scattering in PFN measured around 
Qjv = (0.5,0.5,0.5) Bragg peak. The temperature dependence of the antiferromagnetic Bragg peak is shown in the 

inset. 



